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Abstract 



We present measurements of the excitation function of elliptic flow at midrapidity 
in Au+Au collisions at beam energies from 0.09 to 1.49 GeV per nucleon. For the 
integral flow, we discuss the interplay between collective expansion and spectator 
shadowing for three centrality classes. A complete excitation function of transverse 
momentum dependence of elliptic flow is presented for the first time in this energy 
range, revealing a rapid change with incident energy below 0.4 ^GeV, followed by 
an almost perfect scaling at the higher energies. The equation of state of compressed 
nuclear matter is addressed through comparisons to microscopic transport model 
calculations. 
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The study of collective flow in nucleus-nucleus collisions has been an intense field of re- 
search for the past twenty years [1,2]. At beam energies below several GeV per nucleon, 
the main motivation for studying flow is the extraction of the equation of state (EoS) 
of nuclear matter. This can only be accomplished via comparisons to theoretical trans- 
port models which treat the collision at a microscopic level [3,4]. It is essentially due 
to the competing effects of two-body collisions and mean field dependences that no firm 
conclusion on EoS is established for the moment [4]. 

Elliptic flow at midrapidity (called "squeeze-out" in the early days) is a prominent col- 
lective flow observable that has received great attention in the past. After the pioneering 
measurements at Saturne [5] and Bevalac [6], a wealth of experimental results have been 
obtained at Bevalac and SIS [7,8,9,10,11,12,13,14] as weU as at AGS [15,16], SPS [17,18], 
and RHIC [19]. Elliptic flow is consequently established as a powerful observable in the 
study of relativistic nucleus-nucleus coUisions [20,21]. In what concerns the range of beam 
energies of 0.1-10 AGeV, recent microscopic transport model calculations [4,22,23,24,25] 
have emphasized the importance of elliptic flow for imposing constraints on the mod- 
els, towards the extraction of EoS. In this energy range, the densities, extracted from 
transport model calculations [4], are up to several times normal nuclear matter density. 
Probing EoS at such densities has impfications to astrophysical questions [4]. At beam 
energies below 2 AGeV, no existing set of experimental data can be compared on equal 
footing to the measurements at higher energies. Moreover, recent experimental results 
have demonstrated that directed flow is correlated with transparency [26]. To complete 
the picture, it is important to correlate also elliptic flow with these observations. 

In this paper we present new elliptic flow data at midrapidity for Au-|-Au collisions for 



Corresponding author: GSI, Planckstr. 1, 64291 Darmstadt, Germany; 
Email: A.Andronic@gsi.de; Phone: -|-49 615971 2769; Fax: -|-49 615971 2989. 



2 



the energy range 0.09-1.49 AGeV. The data have been measured with an almost complete 
phase-space coverage using the FOPI detector [27] at SIS, GSI. The reaction products 
are identified according to their charge {Z) in the forward Plastic Wall (PW) at 1.2° < 
9iab < 30° using time-of-flight (ToF) and specific energy loss. In the Central Drift Chamber 
(CDC), covering 34° < 9iab < 145°, the particle identification is on mass (^4), obtained 
using magnetic rigidity and specific energy loss. For the beam energies above 0.4 AGeV, 
measured in a separate nm, the forward drift chamber Helitron, covering the interval 7° < 
9iab < 29° is employed for particle identification on A. We use normalized center-of-mass 
(cm.) transverse momentum (per nucleon) and rapidity p^*^^ = {Pt/A)/{Pp^'/^p)-, V^^^ ~ 
{y/ypY'^': where the subscript P denotes the projectile. Our midrapidity selection is 
|l/(o)| < 0.1. Most of our results are for Z=l particles. For the beam energies above 
0.4 AGeV we also present integral elliptic flow for protons. Identifled by CDC and Helitron, 
pions are excluded from the Z=l sample. For the acceptance of CDC {pf^ > 0.8, at 
midrapidity) ^He is included in the Z—1 sample. For the centrality selection we use the 
charged particle multiplicities, classified into five bins. Ml to M5. We present results for 
the centrality bins M2, M3, and M4, which correspond on average to the geometric impact 
parameter intervals 7.5-9.5 fm, 5.5-7.5 fm, and 2-5.5 fm, respectively. Variations of these 
ranges for different energies, are up to 5% for the energies above 0.4 AGeV and up to 8% 
for the other energies. 

The reaction plane is reconstructed event by event using the transverse momentum 
method [29]. All charged particles detected in an event are used, except a window around 
midrapidity {\y^^^\ <0.3) to improve the resolution. The correction of the extracted V2 
values due to the reconstructed reaction plane fiuctuations is done using the recipe of 
OUitrault [30]. The correction factor is energy- and centrality-dependent. For the M3 and 
M4 centralities this factor is as low as 1.2 for energies around 0.4 AGeV, reaches 1.6 at 
1.49 AGeV and is between 2 and 4.5 for the energies of 0.12 and 0.09 AGeV. For the 
centrality M2, the correction factor is about 1.6 for energies around 0.4 AGeV, reaches 
2.8 at 1.49 AGeV and 4 at 0.12 AGeV. 

Elliptic fiow is quantified by the second order Fourier coefficient V2 — (cos(20)), where (f) 
is the angle with respect to the reaction plane. The systematic errors are for most of our 
energies of the order of 10% for M2 and M3 bins and slightly larger for M4. They increase 
with the beam energy, reaching at 1.49 AGeV about 20%. Except otherwise stated, in 
the following plots, the errors are the systematic ones. They have three sources: i) the 
interpolation over the region 6iab = 30° — 34°, corresponding to p^^^ =0.6-0. 8 at midrapidity, 
which is not covered by FOPI detector, ii) pion contamination in the Z=l sample, as for 
P('^'*<0.2 the pions are rejected on a statistical basis only, iii) particle losses due to multiple 
hits, which occur preferentially in the reaction plane due to the strong directed flow [28]. 

In Fig. 1 we present, for Z—1 particles, the excitation function of elliptic flow integrated 
over pf'^ and of the average values of normalized transverse momentum, (pf'^ ) , for three 
centrality bins. The elliptic flow evolves from a preferential in-plane, rotational-like, emis- 
sion (f2 > 0) to an out-of-plane, or "squeeze-out" (^2 < 0) pattern, as the energy increases. 
As seen in Fig. la, this transition [10,12] depends on centrality. Beyond the transition. 
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Fig. 1. Excitation function for three centrality bins for: a) elliptic flow integrated on momentum; 
b) average normalized transverse momentum. 

for higher energies, the strength of the collective expansion overcomes the rotational-like 
motion [11]. This leads to an increase of out-of-plane emission towards a maximum at 
0.4 AGeV. Beyond 0.4 AGeV, elliptic flow decreases towards a transition to in-plane pref- 
erential emission [15]. In our energy range, the energy dependence of elliptic flow is very 
similar to that of directed flow [1,2,26]. As a consequence of comparable passing times 
and expansion times, in this energy regime, elliptic flow results as an interplay between 
fireball expansion and spectator shadowing. At sufficiently high incident energy, out-of- 
plane emission is expected to grow with: i) the achieved pressure in the fireball (leading 
to larger values, as seen in Fig. lb), which is larger for more central collisions ; ii) 

the degree of shadowing, which is proportional with the impact parameter. This second 
behavior is observed in the data for energies of 0.4 AGeV and above. Fig. la. We wish to 
stress that this phenomenon is quite different from what is known at much higher ener- 
gies, where higher impact parameters lead to increased in-plane emission [20,21,19,18]. In 
this case, corresponding to a much shorter passing time compared to the expansion time, 
shadowing by still incoming participants and by spectators takes place only in the very 
early stages of the collision. The later, dominant, stages of flow generation are determined 
by the almond shape of the isolated fireball [20]. 

In our energy domain, due to increasing transparency in smaller fireballs [26], the growing 
of the out-of-plane emission towards more peripheral collisions will be partially quenched. 
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To elucidate this, it is useful to gather information on generated transverse momenta. 
As seen in Fig. lb, the scaled transverse momenta, (pf^) decrease both as a function of 
centrality and as a function of incident energy. In the low energy end of the data, some 
of the effect as a function of energy is likely due to the decreasing importance, in scaled 
units, of Fermi motion, while at the high energy end, increasing transparency is suggested 
by stopping systematics [26]. In this domain, it is the decrease of {pf'^) which mainly 
determines the decrease of V2 as a function of energy (see next Figure). 




p,(0) 



Fig. 2. Excitation function of differential elliptic ffow for Z=l particles, for M3 centrality. The 
lower and higher energies are shown separately in the left and right panels, respectively. The 
symbols represent the measurements (the errors are statistical) and the lines are parabolic fits. 
In the region not covered by the detector (p^^^^ =0.6-0.8, shades) the symbols are the result of 
the fit. 

More details on the energy dependence of elliptic flow are revealed by a study of the 
transverse momentum dependence of V2, the so-called differential flow. In Fig. 2 we show 
for Z=l particles the differential V2 for all incident energies, for centrality bin M3. In the 
low energy region, differential elliptic flow exhibits a dramatic change as a function of 
beam energy. It evolves from a coexistence of in-plane and out-of-plane features towards 
a monotonic (parabolic) dependence as a function of pf^\ This change may be a result 
of a gradual transition of the collision regime from a mean field to a nucleon-nucleon 
scattering domination. As seen in Fig. 2, the transition from in-plane to out-of-plane 
depends on pf'\ a feature which we have studied earlier in detail [12]. This is markedly 
different than in case of the transition (at about 4 AGeV) from out-of-plane to in-plane, 
which was found to be independent on pt [15]. For energies beyond 0.4 AGeV, our data 
show that elliptic flow approximately scales as a function of pf'\ a feature noted earlier 
for low-pf^ FOPI data [10]. In view of this scaling, it is the decreasing {pt^^) as a function 
of energy (Fig. lb) which explains the decrease of the absolute magnitude of the integral 
elliptic flow for energies above 0.4 AGeV. The scaling of differential V2 is also supported 
by proton differential elliptic flow measured by KaoS collaboration [8] and was noted as 
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well in the case of pions [1]. We interpret it as a consequence of an approximate matching 
between the time scale of the expansion and of the passing time of the spectators, tpass- 
In a simple geometric participant-spectator model, tpass = '^R/ilsVs), where R is the 
radius of the nucleus at rest, Vg is the spectator velocity in cm. and 7^ the corresponding 
Lorentz factor. In the energy range 0.4-1.49 AGeV, tpass decreases from 30 to 16 fm/c, 
implying that overall the expansion gets about two times faster in this energy range. This 
is also supported by the average expansion velocities extracted from particle spectra [7]. 
High-pl'''* particles, for which elliptic flow is larger, are emitted at times even shorter than 
tpass- Originating from the largest density region, they probe the momentum dependent 
interactions of the mean field in a unique way [22]. In a recent paper [13] we have shown 
that detailed time-like diagnosis on the expansion dynamics can be obtained from the 
azimuthal dependence of the collective flow energy. We note that emission times of 13 
fm/c were extracted for high-energy pions at spectator rapidity in Au+Au collisions at 
1 AGeV [31]. 
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Fig. 3. Excitation function of integral elliptic flow for the M3 centrality bin. The FOPI data 
for Z=l particles are compared to measurements available at all fixed target experiments: EOS 
and E895 [15] for protons, E877 [16] and CERES [17] for all charged particles and NA49 [18] for 
pions. The inset shows V2 data for Z <2 particles, analyzed in the directed flow reference frame 
(see text). 

An excitation function of the integral elliptic flow is presented in Fig. 3 for all existing 
measurements at fixed-target experiments. Our data for Z=l particles are compared to 
the data of experiments EOS and E895 [15] for protons, E877 [16] and CERES [17] 
for all charged particles and NA49 [18] for pions. The centrality classes for the various 
experiments are similar to our M3 centrality bin, which corresponds to the range 15%-29% 
of the geometrical cross section. The composite particles included in our Z=l sample are 
important only for energies below 0.8 AGeV, while above, due to their low yield, they 
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do not contribute anymore to enhance the Z=\ flow with respect to proton (see next 
Figure). The increase of elhptic flow beyond 10 AGeV, which is logarithmic in ^/snn up 
to RHIC energies [32], is determined by the gradient-driven expansion of the unhindered 
almond-shape flreball [20]. 

The inset in Fig. 3 shows an excitation function of V2 in the reference frame with the z axis 
(which is usually the beam direction) along the main axis of the directed flow. To allow 
for a direct comparison to the results of Plastic Ball experiment at Bevalac [6] , the results 
are for Z <2 particles, integrated on pf'^ Our data is in a very good agreement with the 
Bevalac results and also agree very well with recent data measured by Indra collaboration 
[14]. The Indra data are for the impact parameter range 4.5-6.5 fm, which is slightly more 
central than the FOPI and Plastic Wall centrality bin. Note that, as the Bevalac data 
were not corrected for the reaction plane resolution, no correction is employed for the 
FOPI and Indra data either. Since all have a quasi-complete coverage of the phase space, 
the correction is similar for the three experiments. It influences the shape of the excitation 
function, in particular at low and high energies. 




Fig. 4. Measured elliptic flow for Z=l (dots) and protons (diamonds) for M3 centrality in 
comparison to BUU model calculations (lines) of: a) Danielewicz et al. [4] and b) Larionov et 
al. [23] for two values of the compressibility, K (in MeV); c) Persram et al. [24] for free and 
in-medium nucleon-nucleon cross sections; d) Gaitanos et al. [25] for static nuclear matter (LDA) 
and non-equilibrium in the treatment of the collision ( ACNM) . 

To assess the implications of our data for constraining the EoS of nuclear matter, we com- 
pare the measurements with calculations using microscopic transport models. In Fig. 4 
we show the comparison of our integral V2 data for M3 centrality with Boltzmann-Ueling- 
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Uhlenbeck (BUU) models. In the upper panels, the sensitivity to EoS, via the compressibil- 
ity is investigated by comparing to calculations of Danielewicz et al. [4] and Larionov 
et al. [23] Note that, performed for 6=5-7 fm, these calculations are for slightly more 
central coUisions than the data. As both sets of calculations are for protons, together with 
our Z—1 data, we include proton flow when available (E6eaTO>0.4 AGeV). As expected, 
the flow is significantly smaller for protons alone compared to Z=l particles only at lower 
energies. There is sensitivity of elliptic flow to EoS, but it is different as a function of 
energy in the two models: larger sensitivity is seen towards higher energies in case of the 
calculations by Danielewicz et al., while the opposite seems true in case of Larionov et 
al. Both models are close to the data, but no consistent agreement between data and 
calculations exists over all the energy range. In particular, both models ovcrprcdict the 
proton data at 0.4 AGeV . As a consequence, no strong constraint on EoS can be derived 
at this stage. Note that the BUU model of Danielewicz et al. was recently shown to favor 
a soft EoS at 0.4 AGeV, in comparison to FOPI data on the azimuthal dependence of 
kinetic energy of hght particles [13]. 

In the lower panels of Fig. 4, the sensitivity of elliptic flow to the treatment of non- 
equilibrium situation of the nucleus-nucleus collision is investigated through the calcu- 
lations of Persram et al. [24] (for all nucleons) and Gaitanos et al. [25] (for Z=l, using 
a coalescence approach to produce composite particles). Both sets of calculations are 
employing a soft EoS, corresponding to ^^=210 and K—2'iQ MeV, respectively. The self- 
consistent treatment of in-medium nucleon-nucleon cross section in the model of Persram 
et al. [24], Fig. 4c, leads to a decrease of the flow magnitude and a better agreement with 
data, in particular for lower energies. The calculations of Gaitanos et al. [25], Fig. 4d, show 
that there is a large effect when including the non-equilibrium momentum distributions 
in the treatment of the collision (asymmetric colliding nuclear matter approximation, 
ACNM), leading to a very good agreement with the data. The equivalent dynamical EoS 
is much softer compared to static nuclear matter under the same conditions (local density 
approximation, LDA). The difference between the two approximations is even larger than 
the span between a soft and a hard EoS has on elhptic flow. 

In general, when comparing elliptic flow data and model calculations, a special care has 
to be devoted to the centrality selection, as the dependence of V2 on impact parameter 
is rather strong (almost linear for the energies above 0.25 AGeV, see Fig. 1). Another 
essential aspect of any comparison between data and model is the particle type selection 
[28,33]. Any conclusion on EoS could be weakened by the inability of the models to 
reproduce the measured yields of composite particles. Other observables, like directed flow 
[28,34] or strange particle yields [35], which may probe the collision in complementary 
ways, need to be simultaneously reproduced as well. Attempts in this direction have 
already been made [4,23,25]. It is unhkely that a narrow bound on EoS can be derived 
based on one single observable. Moreover, as seen in Fig. 4, the differences between various 
models are still large and need to be understood. 

In summary, we have presented new measurements of elliptic flow in Au-|-Au collisions 
for beam energies from 0.09 to 1.49 AGeV. Our data span the onset of the collective 
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expansion, through the interplay of this expansion with the spectator shadowing, which 
we have studied for three centraUty classes. Differential elliptic flow is presented for the 
first time over this broad energy range. Interpreting the observed scaling of elliptic flow 
on transverse momentum for energies above 0.4 AGeV as resulting from an approximate 
matching between the expansion time and the passing time of the spectators, we have 
extracted time scales for particle emission. Our comparison of data and model calculations 
showcase once more the difficulties that need to be solved towards extracting EoS, a 
goal towards which elliptic fiow is a particularly important observable. The precision 
of our V2 data could in principle allow for constraining EoS. It is important that the 
models shall be able to reproduce all the features of V2, including its transverse momentum 
dependence, as well as the corresponding Further improvements and validations of 

current transport models through comparisons with experimental data are needed towards 
a final conclusion on EoS of compressed nuclear matter. 

This work was partly supported by the German BMBF under Contracts No. 06HD953 
and No. RUM-99/010 and by KRF under Grant No. 2002-015-GS0009. 
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